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LATTICE  SYMMETRIES  AND  THERMAL  EXPANSION 


G.  M.  Gallatin 
Perkin-Eiiiter  Coip. 

100  Wooster  Heights  Road 
Danbury,  CT  06810 


ABSTRACT 

It  is  well  known  that  many  of  the  themtotnechanical  properties  of  crystals  are  due  to  anliannonic  effects  in  tlie 
motion  of  the  atoms.  The  anharmonic  terms  in  the  potential  energy  which  lead  to  this  behavior  depend,  in  part,  on 
the  symmetry  of  the  crystal  lattice.  Thus  die  thermomechanical  properties  of  the  crystal,  including  the  coefficient  of 
tliermal  expansion,  depend  on  the  point  group  of  die  lattice,  hi  order  to  study  diis  symmetry  dependence  tlie  path 
integral  representation  of  the  density  matrix  in  the  adiabatic  approximadon  is  used  here  to  evaluate  tlie  thermal 
expectadon  value  of  the  atomic  posidons.  From  this  result  the  approximate  scaling  behavior  of  the  atomic  displace¬ 
ment  with  temperature,  mass,  bond  strength  (harmonic  coupling),  and  anharmonic  coupling  is  determined.  The  path 
integral  calculation  is  fully  quantum  mechanical  and  can  be  used,  in  principle,  to  derive  all  tlie  diermomechanical 
properties  of  the  crystal  within  die  adiabadc  approximation. 


1.  INTRODUCTION 

Tlie  standard  approach  for  deriving  many  of  the  properties  of  crystals  is  tlie  harmonic  approximation  which  in  quan¬ 
tum  language  becomes  the  phonon  representation  of  lattice  vibrations  [1].  Unfortunately,  tliis  approximation  does 
not  include  many  of  the  known  ihennomechanical  properties  of  crystals  [1].  For  example,  tlie  size  of  a  rigorously 
harmonic  crystal  is  independent  of  the  temperature  and  thus  in  this  approximation  the  coefficient  of  tliermal  expan¬ 
sion  of  all  crystals  is  zero.  There  have  been  many  studies  of  the  effects  of  anliannonicity  on  crystal  properties 
[l]|2].  They  usually  proceed  by  using  the  phonon  (hitmionic)  representation  of  lattice  vibrations  as  a  basis  for  die 
perturbative  calculation  of  phonon-phonon  (anharmonic)  interaction  effects.  Here  we  use  the  patit  integral  approach 
developed  by  Feynman  for  calculating  the  quantum  density  matrix  [3].  It  provides  a  compact  and  useful  representa¬ 
tion  of  quantum  mechanical  systems  at  finite  temperature.  Using  the  density  matrix  we  can  calculate,  in  principle,  all 
the  lliemioniechanical  properties  of  die  crystal.  Note  that  since  die  padi  integral  uses  die  entire  normal  mode  spec- 
trum  of  the  lattice,  it  includes  both  the  distortion  of  unit  cells  as  well  as  the  motion  of  the  unit  cells  with  respect  to 
each  other. 


2.  PATH  INTEGRAL  FOR  THE  ANHARMONIC  OSCILLATOR 

In  this  section  we  present  the  path  integral  result  for  a  single  anharmonic  oscillator.  This  will  be  generalized  below 
to  cover  the  coupled  oscillators  (atoms)  wltich  describe  a  crystal  lattice. 

As  shown  in  [3]  the  quantum  mechanical  density  matrix  for  a  single  particle  of  mass  M  interacting  with  a  potential 
V(.v,  s)  at  temperature  T  is  given  by 


-  1  - 


(I) 


The  notation  is  as  follows:  |3  =  l /kBT  where  kB  is  the  Boltzmann  constant,  7?  is  Ploncks  constant,  s  is  a  dummy  inte¬ 
gration  variable  which  parametrizes  the  "path"  v(r)  and  acts  as  a  "time”  variable,  and  &v  indicates  integration  over  all 
possible  paths  x(s)  that  satisfy  ,v(0)  =  x,  and  ,t(p7I)  =  ,rr 

All  thermodynamic  quantities  can  be  calculated  from  the  density  matrix.  For  example,  the  partition  function,  Z,  is 
given  by  the  trace  of  the  density  matrix 

Z  =  rr[p]  =  J  dx  p(.t,  .t)  (2 ) 

fhe  potential  C(.v,  s )  can  be  expanded  in  a  power  series 

V(x,  s)  =fis)x  +  —■kx2  +  C3 x3  +  C4 x*  +  ...  *fis\x  +  1-kx  +  U{x)  (3) 

L,  ^ 

where  fis)  represents  an  external  driving  force,  k  is  the  harmonic  oscillator  force  constant  and  die  coefficients 
C3,  Cv  ...  are  the  cubic,  quartic  and  higher  order  anharmonic  coupling  constants,  respectively.  Utx)  is  called  the 
anharmonic  potential. 

Using  a  standard  trick  [3]  we  can  write 

P(.ta,  ■',)  =  expf"-  ds  U(- ITS/Sfa))!  | exp|^-  -jr  j ds  ( y  (dx)1  +  y  .tz  J 

■eXp[-{|t']P»W  <4) 

in  an  obvious  notation.  Note  that  p„  is  a  functional  of  the  external  force  fix).  Expanding  the  first  exponential  in  a 
power  series  yields  the  perturbation  expansion  of  the  density  matrix.  The  functional  derivatives  8/5ftj)  can  be  per¬ 
formed  once  U  is  written  as  the  anharmonic  series  in  (3).  In  many  cases  fis)  =  0.  Thus  after  performing  Hie  function¬ 
al  derivatives  to  obtain  the  perturbation  to  a  given  order,  fis)  should  be  set  to  zero.  Clearly,  in  any  case  where  tlte 
external  force  is  not  zero  we  set  fis)  equal  to  die  given  external  force  after  taking  the  functional  derivatives. 

p0  can  be  evaluated  exactly  using  the  tecltnique  of  stationary  phase  [3].  Tlie  result  is 

where  £„[/]  is  the  "action"  for  the  "classical"  trajectory  with  co=  Vk/Af.  The  reader  is  referred  to  Ref.  [3]  where  50 
is  given  explicitly. 

Since  we  wish  to  study  the  thermal  expectation  values  of  position,  <x">  with  n  =  1,  2,  ...  we  compute  the  generat¬ 
ing  function  for  these  expectation  values 
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Z(0) 
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<  „p(  -m) , = mau&l . 
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From  the  definition  of  Z(r |)  we  have 


dr  p(.r,  x) 


n 

<X  > 


[<-a, ran)] 


'T1  a  0 


Z(  0) 


Performing  the  traces  in  (6)  yields 


(7) 


Z(r|)  = 


r  1  f 

rrl 

HJ' 

U  exp 

ia  +  jL  +  i 

2d  4d  2 


J  dsds'fts)  G(s,  s') /(s') 
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where 


(8) 


_  Afto  cosh(p7ioj)  -  1 
7i  sinh(P7ko) 

_  Afto  (/  +  /)  (exp(p7ico)  -  1) 

7i  sinh(pnm) 

/=•——  f  ds  exp(  -  Kras)  f{s) 

2M(0  J 
0 

1  =  2^to  1 ds  exp(  “  T" ~ 

0 

I 

G(s,  s')  =  -  —77 —  exp(  -  7ico)s  -  s'l) 

2A/co 

Using  the  above  results  we  can  evaluate  <  x  >  (0  lowest  older 


(9) 


<  x  >  =*  -  C - — —  coth(P7ko/2) 

2 MV 


(10) 


•4 


This  is  the  standard  result  [2][3]  and  shows  how  the  average  position  of  the  particle  (to  lowest  order*  depends  on  the  D 

temperature.  Note  that  <x>  does  not  go  to  zero  at  T  =  0  due  to  (he  quantum  zero  point  fluctuations  of  die  panicle  O 

position.  - 
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3.  THE  ANHARMONIC  CRYSTAL 

The  atoms  io  a  crystal  can  be  treated  as  coupled  anharmonic  oscillators.  In  the  adiabatic  approximation  f  1 1  the  elec¬ 
tron  degrees  of  freedom  are  replaced  by  an  effective  potential  energy  V  which  depends  on  tlie  positions  of  llie  atoms. 
The  classical  zero  temperature  equilibrium  position  of  each  atom  is  detenniued  by 


3V  =  0  (11) 

For  a  crystal  the  equilibrium  positions  form  a  lattice  of  points  that  can  be  labelled  by  a  unit  cell  index, 
m  =  mv  mr  m3  and  a  basis  index  p  (i  |[2j.  For  a  primitive  lattice  with  only  one  atom  per  unit  cell,  p  =  1.  whereas  for 
a  non-primitive  lattice  with  S  atoms  per  unit  cell  p  =  1,...,S.  At  any  finite  temperature  the  atoms  vibrate  around 
their  equilibrium  positions.  Generally  these  vibrations  are  not  too  large  and  we  can  expand  V  around  the  equihbiium 
positions,  R"  ,  where  a  labels  the  xyz  components  of  the  position  vectors.  For  ease  of  notation  we  can  collectively 
denote  the  indices  m  p  a  by  the  index  A.  Letting  xA=RA  +  qA,  the  expansion  of  V  about  the  equilibrium  positions,  RA 
takes  the  form 


V  qA  +  1  X  CA*  Wb  +  X  CASC  Wc  +  •••  <12> 

A  AB  ABC 

As  before  the  anharmonic  terms  in  tlie  above  expansion  will  be  denoted  by  the  anharmonic  potential  U,  llie  expan¬ 
sion  coefficients,  C,  are  the  harmonic  and  anharmonic  coupling  constants  and  the  functions  fA(s)  are  external  forces. 

For  the  crystal  in  the  adiabatic  approximation  the  argument  of  the  exponential  in  the  path  integral  for  p0  is  -I/7i 
times  the  integral  with  respect  to  s  of 


J  X  Ma{3,^2  +  }£Cab  Wb  +  XfA(s)  qA  ,13» 

A  AB  A 

where  Ma  =  is  the  mass  of  tlie  p  r/i  type  of  atom  in  tlie  basis.  (13)  can  be  written  as  tlie  sum  of  uncoupled  normal 
mode  harmonic  oscillators  using  the  transformation  11] 


Os) = X  s) 

ka 


(14) 


with 


r  a)  =  — = 


eMa  (k)  exi 


(15) 


where  e  is  the  polarization  vector  of  the  normal  modes  labelled  by  a,  N  the  number  of  atoms  in  the  crystal,  and  k 
tlie  spatial  frequency.  In  obtaining  (15)  we  have  used  tlie  standard  Bom-von  Kannan  (periodic)  boundary  conditions 

UJ. 


Substituting  (14)  into  (13)  gives 
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i  I  d£a(l  s)  |2  +  y  0) \Ck)  |  %a(l  s)  f  +  gfk,  s)%a(t  *)J  (16) 

where  «„<*)  is  the  fundamental  frequency  of  the  £  a  normal  mode  and  ga(k,  s)  =  \  rAa(k)f^s). 

With  this  factorization  p  and  Z(r|a(fc))  take  the  same  form  as  in  (4)  and  (8)  but  with  a  sum  in  the  exponent  over  all 
the  k  a  values.  As  in  the  single  panicle  case  the  expectation  value  of  die  position,  qA,  can  be  evaluated  older  by 
order  in  the  anhaimonic  potential  by  taking  the  appropriate  derivatives.  Tlie  lowest  order  contribution  from  the  cubic 
coupling  term  takes  the  form 


I 


1 

8(toa(0))2 


(17) 


where  C  is  given  by 


c,<*-)  -  2  CA»c  * «•  ]  (II SI . 

ABC 

with  c.c.  being  the  complex  conjugate  of  the  proceeding  term. 

In  the  above  result  only  optical  phonons  (  coa(0)*0  )  contribute.  Therefore  this  term  is  non-zero  only  for  non¬ 
primitive  lattices.  For  primitive  lattices,  higher  order  terms  in  the  expansion  must  be  evaluated.  This  can  be  done 
straightforwardly  by  computing  the  higher  order  functional  derivatives  in  the  perturbation  expansion  of  Z.  Note  that 
(17)  has  the  same  form  as  (10)  as  would  be  expected.  Again  due  to  the  quantum  zero  point  fluctuations  <qA>  does 
not  vanish  at  T  =  0. 

Detailed  predictions  from  (17)  and  higher  order  terms  require  the  input  of  lattice  symmetries  and  specific  coupling 
constant  strengths  as  discussed  in  the  next  section.  But,  we  can  use  (17)  to  determine  the  approximate  settling  depen¬ 
dence  of  <qA>  with  temperature,  mass,  bond  strength  (harmonic  coupling)  and  anhanuonic  coupling.  Since  we 
want  an  approximate  scaling  law  we  will  ignore  the  £  a  dependence  of  to  and  M.  Note  that  with  tltis  assumption  we 
are  losing  tlie  distinction  between  the  mass  of  single  atoms  and  tlie  reduced  mass.  For  high  temperatures,  i.e..  temp¬ 
eratures  with  k^T »Tmma  where  co^  is  the  maximum  optical  phonon  frequency,  coth( p/ico/2 )  -  l/|)7ico.  Taking 

at  -  V  KIM  where  K  is  die  harmonic  force  constant  or  bond  strength  and  including  the  M  dependence  of  T,  we  find 


<qA>~ 


CT 

K2M 


(19) 


This  result  is  very  similar  to  other  scaling  laws  ami  shows  a  dependence  on  mass,  coupling  and  temperature  which  is 
in  accord  with  expectation,  and  matches  that  predicted  by  equation  (10). 
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4.  LATTICE  SYMMETRIES 


Uniform  translation  and  rotation  of  the  crystal  as  well  as  the  symmetry  of  the  lattice  will  impose  restrictions  on  the 
coupling  constants.  These  requirements  have  been  worker!  out  in  detail  for  many  simple  crystals.  Here  we  present 
only  a  few  of  these  results.  For  more  details,  see  Refs.  [l][2]. 

Since  the  coupling  constant  are  given  by  derivatives  of  V  evaluated  at  the  classical  zero  temperature  equilibrium 
positions  and  derivatives  commute,  the  coefficients  must  be  invariant  under  interchange  of  their  indices,  i.e., 


r  =C  =  C 

ABC...  BAC...  CBA. 


(20) 


Uniform  translation  of  the  crystal,  which  is  equivalent  to  the  conservation  total  linear  momentum,  yields  die  con¬ 
straint 


'2" 

rtV 

Uniform  rotation  of  the  crystal,  which  is  equivalent  to  the  conservation  of  total  angular  momentum,  yields  con¬ 
straints  which  relate  the  n  and  n  +  1  order  constants.  The  relation  between  die  quadratic  and  cubic  coupling  con¬ 
stants  is 


,  mnl 


\»avhkc’ 1  Xc 


5  +C' 


.6 

JUJV;  be 


(K 

is  symmetric  under  the  interchange  of  c  and  c. 

The  symmetries  that  depend  specifically  on  the  lattice  yield  constraints  that  take  the  form 


(22) 


C 


ga...V* 


-2«>. 


M  c...vd 


cd 


(23) 


wliere  tlie  matrix  O  is  an  element  of  the  lattice  point  group.  Tlius  £1  maps  the  positions  labelled  by  nip . .  ,nv  to  (hose 
labelled  by  nip...nV. 

These  symmetry  constraints  determine  which  elements  of  the  coupling  constants  are  zero  and  what  the  relationship 
is  between  the  remaining  non-zero  elements.  The  absolute  value  of  die  non-zero  elements  Ls  determined,  of  course, 
by  the  strength  of  the  interatomic  bonds  which  are  in  many  cases  only  approximately  known.  But,  even  in  the 
absence  of  accurate  values  for  these  bond  strengths  a  range  of  possible  expansion  behaviors  can  be  determined  by 
using  the  symmetry  constraints  and  varying  the  absolute  values  over  reasonable  ranges.  Furtlier.  tlie  effective  poten¬ 
tials  developed,  for  example,  in  Refs  [4)15],  can  be  used  to  directly  calculate  the  coupling  constants  and  thus  evalu¬ 
ate  tlie  thermal  position  dependence  of  the  atoms. 
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5.  SUMMARY 

The  path  integral  technique  for  evaluating  the  quantum  density  matrix  has  been  used  to  detennine  the  form  of  ilte 
thermal  expectation  value  of  the  position  of  an  atom  in  a  crystal  lattice.  From  this  result  we  have  predicted  llte  low¬ 
est  order  scaling  behavior  of  the  atomic  position  with  mass,  temperature,  bond  strength  and  anhamtonic  coupling. 
This  scaling  law  is  in  accord  with  expectation.  Further,  the  result  explicitly  contains  the  anliarmonic  coupling  con¬ 
stants  which  depend  on  the  symmetry  of  die  lattice.  Thus,  using  this  formalism  it  is  possible  to  directly  calculate  the 
effects  of  lattice  symmetry  ou  die  thermal  position  dependence  of  the  atoms  in  the  lattice.  From  these  results  it  is 
possible  to  determine  a  range  of  thermal  expansion  behaviors  for  each  type  of  lattice.  The  results  of  otlier  tlieoreti- 
cal  analyses,  such  as  those  in  [4]  and  15]  can  be  substituted  into  the  above  fonnulation  to  obtain  predictions  of  all  the 
thennomechanical  properties  of  each  type  of  lattice. 
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